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Abstract:In order to extend the applications GNSS-R/IR (GNSS-Reflectometry 

/Interferometric Reflectometry ) remote sensing technique to soil freeze/thaw process detection, 
soil(frozen/thawn) mixing permittivity models are employed to calculate the soil permittivity. 
Bistatic full-polarizationcoherentscattering model and random roughness surface scattering model 
are used to calculate the coherent and non-coherent scattering, which result in the variations of GPS 
multipath observables and DDM (Delay Doppler Map) waveforms, respectively. When the soil 
freeze/thaw process occurs, theoretical simulations indicate that soil freeze/thaw process induce the 
abrupt permittivity changes and affect the obvious variations of GNSS-R/IR signals.In this way, 
theoreticalfundamentals for soil freeze/thaw process detections are presented. 
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摘要 : 
GNSS-R/IR (GNSS-Reflectometry/Interferometric Reflectometry) 遥 感 是 利用 导航 卫星 反射 信号 或 
直射 信号 与 反射 信号 的 相干 信号 对 地 观测 的 新 兴 遥 感 方式 ， 在 国内 外 引起 广泛 关注 。 为 将 GNSS- 
R/IR 技术 的 应 用 领域 拓展 到 地 表 冻 融 状 态 的 监测 上 ， 本 文 利 用 冻 融 土 混合 介质 介 电 常 数 模型 计 
算 土 壤 介 电 和 常数， 采用 双 站 全 极 化 相干 反射 率 模型 和 随机 粗糙 面 散射 模型 ， 分别 计算 了 经 冻 融 
土 反 射 的 GPS 相干 反射 量 的 镜像 反射 率 , 以 及 GPS 非 相 干 反射 分 量 的 漫 散 射 特性 ; 并 模拟 分 析 了 
冻 融 转换 时 ，GPS 多 路 径 信息 CGNSS-IRO 以 及 包含 漫 散 射 信号 的 延迟 多 普 勒 图 (GNSS-R) 的 变化 
特征 。 理论 研究 表明 冻 融 转换 过 程 中 , 地 表 介 电 常 数 的 变化 导致 GPS 多 路 径 信息 和 延迟 多 普 勒 图 
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的 明显 变化 。 本 文 从 散射 机 理 上 揭示 了 利 
论 依据 。 


| GNSS-R 和 GNSS-IR 遥感 进行 地 表 冻 融 特性 监测 的 理 


关键 词 :GNSS-R/IR， 冻 融 地 表 ， 延 迟 多 普 勒 图 ， 多 路 径 ， 双 站 全 极 化 
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季节 性 冻 土 和 永久 性 冻 土 约 占 地 球 低 成 本 、 小 功 耗 、 轻 体积 以 及 高 时 空 分 辩 


陆地 总 面积 的 33%， 主 要 分 布 在 高 纬度 和 
高 海拔 地 区 。 陆 地 表层 土壤 冻 融 状态 转换 
随 季节 每 年 重复 发 生 ， 与 人 类 生活 环境 密 
切 相 关 ， 由 于 土壤 中 水 的 相 态 变化 ， 这 一 
过 程 强烈 地 影响 着 地 表 辐 射 能 量 的 转换 过 
程 和 地 表 径 流 的 强度 ， 是 地 表 能 量 平衡 和 
水 分 平衡 的 重要 影响 因素 ,也 是 气候 变化 的 
灵敏 指示 器 ,因此 有 效 监 测 地 表 冻 融 状 态 的 
时 空 分 布 及 其 相关 物理 参数 是 冰冻 圈 、 地 
学 、 水 文学 等 研究 领域 的 重要 课题 。 

可 见 光 和 热 红 外 遥感 受 天 气 条 件 限 
制 ， 微 波 遥 感 却 可 以 全 天 时 全 天 候 的 观 
测 。 主 /被 动 微波 遥感 (雷达 /辐射 计 ) 是 地 


表 冻 融 状 态 监 测 的 有 效 手 段 之 一 ， 星 载 观 
测 将 空前 提高 地 表 冻 融 监测 的 空间 分 辩 


率 ， 但 其 时 间 分 辨 率 (每 3 天 全 球 重复 覆 
盖 ) 与 实际 监测 的 科学 需求 之 间 存在 一 
的 差距 。 
近年 来 ， 利 用 导航 卫星 的 反射 信号 
(GNSS-R) 或 直射 信号 和 反射 信号 的 相干 信 
号 (GNSS-IR) 进 行 陆 面 参数 遥感 研究 成 为 一 


种 新 兴 的 对 地 观测 方式 ，GNSS-R/IR 显著 的 
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等 诸多 优点 可 以 成 为 现 有 地 面 站 点 观测 
和 机 载 / 星 载 遥感 监测 的 有 效 有 益 补 充 ”。 
GNSS-R 的 工作 模式 分 为 散射 计 模式 和 
高 度 计 模式 ， 针 对 GNSS-R 散射 计 模 式 ， 需 
要 研制 专门 GNSS 反射 信号 接收 机 对 陆 面 
参数 进行 遥感 监测 ， 在 该 种 接收 机 的 基础 
上 利用 其 相关 功率 值 遥感 地 物 参数 。 目 前 
陆地 表面 研究 领域 涉及 土壤 水 分 、 植 被 售 
水 量 和 生物 量 以 及 积 雪 厚度 研究 。 
GNSS-IR 遥感 不 需要 研制 专门 的 接收 
机 ， 可 以 直接 使 用 地 球 测绘 或 者 地 球 物 理 
的 接收 机 对 地 表土 壤 水 分 、 植 被 参数 、 积 


雪 厚 度 和 潮 位 变化 进行 遥感 监测 4 |, H 
空间 分 辨 率 约 为 lkm: 介 于 传统 站 点 式 传 


感 器 («1m^ 和 星 载 观测 (2100km^ 之 
间 。 利 用 从 多 路 径 观测 数据 中 的 有 效 反射 
计 高 度 ， 幅 度 和 相位 信息 可 以 提取 相关 地 
物 参 数 。 美 国 PBO H20 研究 团队 将 上 述 环 
境 遥 感 产品 每 隔 12 小 时 公布 一 次 中 ， 其 土 
壤 水 分 反 演 结果 已 被 用 做 SMAP (Soil 

Moisture Active Passive) 卫星 的 校 验 中 四 。 


目前 利用 GNSS-R/IR 进行 地 表 冻 融 状 
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态 监 测 的 研究 相对 较 少 ， 
RIR 对 地 表 冻 融 状态 监测 存在 着 可 行 性 
是 一 种 新 兴 的 监测 手段 。 研 究 者 已 经 利用 
相关 理论 和 IGS (International GNSS 
Service) 站 点 数据 分 析 中 进行 了 一 些 探索 
性 研究 ""10。 本 文 详细 阐述 了 利用 GNSS- 
RIR 遥感 方式 进行 地 表 冻 融 监 测 的 理论 机 
理 。 

针对 星 载 观测 ，GNSS 卫星 群 和 专门 
的 GNSS-R 反射 信号 接收 机 之 间 形 成 一 种 
双 站 雷达 的 工作 模式 ， 接 收 机 得 到 的 DDM 
波形 是 用 于 反 演 地 球 物理 参数 的 基本 观测 
量 。 开 发 建立 相应 的 仿真 软件 是 开展 卫星 
观测 、 数 据 仿真 、 实 验 设 计 ， 内 插 外 推 数 
据 和 反 演 等 工作 的 重要 关键 机 理工 具 
n223。 目 前， 国内 外 的 研究 团队 重点 针对 
海 况 参数 展开 仿真 分 机， 已 经 发 展 建立 


CYGNSS E2ES (cyclone global navigation 


但 利用 GNSS- 


satellite systems mission end to end simulator 
(E2ES)) 软件 04，P2EPS (the PAU/PARIS 
end-to-end performance simulator) 仿 


真 软件 03 和 GREEPS (an GNSS-R end-to- 


end performance simulator) 仿真 软件 
H6] SAVERS (a simulator of GNSS 
reflections from bare and vegetated 
Soils) 仿真 软件 则 针对 陆 面 土壤 和 植被 特 
性 研究 而 建立 WT。 近期 ， 研 究 人 员 在 此 Z- 
V 模型 基础 上 MW， 发 展 建立 了 地 表 冻 融 转 
换 DDM 仿真 分 析 模 型 (31。 

大 多 数 的 前 向 GPS 多 路 径 模型 关注 码 
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调制 ， 对 于 反射 能 量 、 相 位 和 延迟 采用 任 
意 值 ， 或 者 基于 给 定 几 何 计算 反射 延迟 ， 

亦 或 根据 经 验 值 定义 反射 能 量 。 基 于 物理 
模型 的 前 向 GPS 多 路 径 模型 ， 最 早 的 是 

Zavorotny 等 人 的 模型 tI， 该 模型 完全 考虑 
了 测绘 接收 机 和 天 线 的 直射 信号 和 反射 信 
号 的 极 化 特性 ， 采 用 分 层 的 办 法 计算 反射 
言 号 ， 并 利用 小 波 扰 模 型 计算 反射 率 ， 将 
介 电 常数 转换 为 反射 系数 ， 最 后 把 反射 系 
数 和 指定 的 天 线 辐射 样 图 耦合 得 到 SNRmpi 


(multipath Signal Noise Ratio) 干涉 图 。 基 


于 Zavorotny 等 人 的 物理 模型 ，Nievinski 
和 Larson 给 出 了 一 个 全 极 化 的 前 向 GPS 
多 路 径 模型 20， 该 模型 考虑 了 GPS 广播 信 
号 的 右 旋 和 左旋 圆 极 化 信息 、 天 线 朝向 、 
地 表 响 应 、 伪 码 测 距 和 噪声 能 量 。 本 文 在 
此 模型 基础 上 ， 发 展 建立 了 地 表 冻 融 状 态 
转换 模型 。 

本 文 将 已 有 的 DDM 仿真 分 析 模型 和 
前 向 GPS 多 路 径 模型 拓展 到 适宜 进行 地 表 
冻 融 状态 监测 的 理论 模型 ， 并 以 此 模型 为 
工具 ， 模 拟 分 析 了 地 表 冻 融 状态 转换 对 
DDM 波形 和 多 路 径 数据 的 影响 ， 从 散射 机 
理 上 揭示 了 利用 GNSS-R/IR 遥感 进行 地 表 
冻 融 特性 监测 的 理论 依据 ， 为 后 续 利 用 
GNSS-R/IR 方法 进行 地 表 冻 融 监测 提供 理 
论 依据 。 


l. 地 表 冻 融 转 换 时 的 电磁 散射 特性 变化 
物质 的 复 介 电 常 数 是 其 固有 属性 , 实 部 
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反应 了 不 同 介质 表面 发 生 的 波 的 折射 和 反 
射 现 象 ; 虚 部 则 表示 入 射电 磁 波 在 介质 中 的 
衰减 (吸收 和 转化 ) 程度 。 本 部 分 利用 混合 
介 电 常数 模型 模拟 分 析 了 地 表 冻 融 转换 时 ， 
介 电 常数 以 及 其 对 应 的 各 种 极 化 条 件 下 的 
反射 信号 的 相干 和 非 相 干部 分 的 散射 特性 
变化 情况 。 
1.1 冻 融 转换 导致 介 电 常数 差异 

基于 电磁 波 传输 理论 ,冻结 土壤 可 看 做 
是 由 空气 、 固 体 颗粒 、 自 由 水 、 束 缚 水 和 冰 
五 种 物质 组 成 的 介 电 混 合体 ， 其 复 介 电 常 数 
理论 上 受 入 射电 磁 波 的 频率 、 温 度 、 土 壤 含 
水 量 等 因素 影响 。 最 终 介 电 常数 可 看 做 是 各 
组 分 介 电 常数 之 和 ， 如 公式 1 所 示 。 
Ea 一 > Ver a» 


其 中 ，% 是 形状 常数 因子 ，Vi 是 第 i 个 组 分 


的 体积 百 分 含 量 。 


本 文 土壤 介 电 常 数 的 实 部 
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和 虚 部 可 通过 普遍 采用 的 半 经 验 模型 一 一 
Dobson 介 电 常数 模型 计算 得 到 它 约 。 图 1 
给 出 利用 该 模型 模拟 的 不 同 土壤 水 分 含量 
下 ，GPSL1 载波 频率 ， 土 壤 介 电 常 数 实 部 
和 虚 部 随 土壤 温度 的 变化 情况 。 从 中 可 以 
看 出 当 土壤 从 冻结 状态 转换 到 融化 状态 过 
程 中 ， 在 各 种 土壤 水 分 含量 下 ， 介 电 常 数 
的 实 部 和 虚 部 均 变 化 明显 ， 在 土壤 温度 大 
于 0"C 时 ， 土 壤 水 分 含量 对 介 电 常数 影响 
( 实 部 ， 虚 部 ) 明显 ， 在 土壤 温度 小 于 0°C 
时 ， 由 于 冰 的 存在 ， 土 壤 水 分 对 实 部 有 影 


响 但 不 明显 ， 对 于 虚 部 的 影响 完全 可 以 和 


略 。 
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Soil Temperature (C) Soil Temperature (C) 
1 GPSL1 载波 频率 时 ， 在 不 同 土壤 水 分 含量 下 ， 介 电 常 数 实 部 (a) 和 虚 部 (b) 随 土壤 温度 的 变化 情 
况 


Fig.1 Dielectric constants (real part (a) and imaginary part (b)) versus soil temperature of GPS L1 carrier 


frequency at different soil moisture. 


对 于 不 同 的 GPS 载波 频率 ， 土 壤 的 介 
电 常数 如 表 1 所 示 , 由 表 1 可 知 在 冻 融 转换 
时 ，GPS 载波 频率 ( 均 为 低频 工 波段 ) 对 介 


电 常数 的 影响 可 以 忽略 ,因此 在 后 续 模拟 分 
析 中 只 采用 L1 载波 频率 时 的 介 电 常数 开展 
分 析 。 


表 1 GPSL1， L2 ALS 载波 频率 下 冻 Cl °C) wR (1 °C) 土壤 的 介 电 常 数 


Tab.1 Dielectric constants of GPS L1, L2 and L5 frequency band as soil change from frozen (-1 °C) to 


thawn state (1 °C) 


Ll 载波 频率 L2 载波 频率 L5 载波 频率 
Ts (C) | 实 部 虚 部 实 部 虚 部 实 部 虚 部 
-1 8.66769 1.03403 8.70793 1.16876 8.71293 1.19762 
+1 21.10970 3.54415 21.27701 3.49827 21.29776 3.51016 
1. 2 冻 融 转换 时 圆 极 化 反射 率 变 化 情况 
在 考虑 相位 的 情况 下 , 利用 极 化 合成 的 方法 ， (Ir. í +f) 1 . 
RR 二 Orr = 4 -= Re(rr, ) 


计算 得 到 RR 和 LR 极 化 下 的 反射 率 P529 
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O 为 镜像 双 站 雷达 散射 截面 , 下 角 标 R. 


L 分 别 代 表 RHCP (Right Hand Circular 


BRCS(dB) 


Polarization ) 和 LHCP (Left hand circular 
polarization) 极 化 状态 。7 为 镜像 反射 系数 ， 
TER v M h 分 别 代表 垂直 极 化 和 水 平 极 化 。 


Ele (°) 


图 2. 地 表 冻 (-1°C) 融 (1°C) Bj, [BIÉE(VR/HR) (a) / 圆 (RR/LR) (b) 极 化 BRCS 变化 


Fig.2 Circular-linearly(VR/HR)(a) and circularly (RR/LR)(b) BRCS as soil change from frozen to thawn state 


2 给 出 了 当 土 壤 温 度 分 别 为 -1%C OR 
结 ) TU lC (融化 ) 时 ， 圆 极 化 和 圆 线 极 化 
在 同 极 化 和 交叉 极 化 下 的 BRCS (Bistatic 
Radar Cross Section) 变化 情况 。 由 图 可 知 ， 
在 布鲁斯 诺 角 附近 ，VR 极 化 存在 BRCS 最 
低 点 ， VR 极 化 在 小 仰角 ， 即 1°C 布 鲁 斯 诺 
角 之 前 ，-1°C 的 BRCS 高 于 1°C， 而 在 -1°C 


布鲁斯 诺 角 之 后 ，1°C 时 的 BRCS 高 于 -1°C 
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的 BRCS; 在 图 2 中 给 出 的 仰角 范围 内 , -1°C 
时 的 BRCS 高 于 1"C 的 BRCS (两 条 BRCS 
交叉 线 之 间 ) ， 在 BRCS 交叉 线 之 后 和 第 二 
个 布鲁斯 诺 角 之 前 , IC 的 BRCS 高 于 -1°C，; 

RR 极 化 随 着 仰角 的 增加 ，BRCS 增加 ，1°C 
的 BRCS 高 于 -1"C 的 BRCS; LR 极 化 的 
BRCS 随 着 仰角 增加 而 降低 ，-1*C 的 BRCS 


高 于 1°C 的 BRCS。 
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1.3 冻 融 转换 时 BRCS 变化 


实际 地 表 为 随机 粗糙 表面 ， 入射 的 导航 
卫星 信号 经 地 表 散 射 后 , 散射 能 量 中 的 非 相 
干部 分 不 可 忽视 。 而 常用 的 随机 粗糙 地 表 的 
面 散射 模型 有 KA (Kirchhoff Approach) 模 
型 、 SPM (Small Perturbation Method) 模型 、 
IEM (Integrated Equation Model) 模型 以 及 
AIEM (Advanced Integrated Equation Model) 
模型 。 AIEM 模型 是 在 IEM 模型 的 基础 上 的 
改进 模型 ， 形 式 更 为 复杂 ， 但 是 模型 的 精度 
得 到 了 提高 ，AIEM 模型 的 基本 形式 可 以 概 
括 为 三 项 之 和 : 基 尔 霍 夫 项 a* ， 基 尔 霍 夫 


VV FTBRCS (dB) 


E 


HH FT BRCS (dB) 
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补偿 项 os 和 二 者 的 交叉 项 o% ， 如 公式 4 
所 示 。 公 式 中 o^ 为 双 站 雷达 散射 截面 ，q、 
p 分 别 为 接收 和 发 射 时 的 极 化 状态 27。 

dx Sues (4) 

图 3 是 利用 极 化 合成 方法 计算 得 到 的 地 
表 冻 融 转换 时 各 种 极 化 BRCS 的 变化 情况 。 
在 地 表 从 冻结 到 融化 转换 的 过 程 中 ，BRCS 
在 各 种 极 化 和 角度 时 ， 变 化 差异 明显 ， 即 观 
测 几何 中 的 天 项 角 和 方位 角 以 及 观测 的 各 
种 极 化 都 会 导致 冻 融 转换 时 BRCS 的 变化 。 
如 何 从 角度 和 极 化 信息 中 有 效 提取 冻 融 转 
换 敏感 参数 是 后 续 研 究 重点 内 容 之 一 。 


VH FT BRCS (dB) 


-0.5 -0.5 
[7] e^ [7] 
全 E = 
O 0 o 0 Ng 
[z 2 = 
[7] o [7] 
0.5 0.5: 
1 1 ; ; 
-1 0 1 -1 0 1 -1 0 1 
| (a) 0) . (c) 
sind Coso, sind cosó. sind Coso, 
RR FT BRCS (dB) RL FT BRCS (dB) 


s 
LI 


sind sing 
sind sing 


(d) 


sind. COSÓ, 


图 3. 冻 (-1°C) 融 (1°C) 转 换 时 各 种 极 化 ( 线 极 化 、 圆 极 化 ) 下 BRCS 变化 差异 


(e) 


sin@_cosé 
S s 


Fig.3 BRCS differences as the soil temperature changes from frozen state (-1°C) to thawn state (1°C) at 


linear polarization (VV, HH and VH) and circular polarization (RR and RL) 
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2 . 冻 融 转化 时 GNSS-R/IR 接收 机 信号 变化 


2.1 多 路 径 信息 变化 


当地 物 参数 , 如 积 雪 厚度 、 植 被 含水 量 、 
土壤 水 分 变化 时 , 会 引起 地 球 测绘 或 者 地 球 
物理 的 GPS 接收 机 的 多 路 径 信 息 的 变化 
(SNR, MAE) 。 因 此 ， 可 以 利用 直 
射 信 号 和 反射 信号 的 相干 多 路 径 数据 对 地 
物 参数 进行 反 演 。 
Nievinski 和 Larson 发 展 建立 了 可 以 同 
时 考虑 GPS 信号 极 化 、 天 线 和 地 表 响 应 的 
全 极 化 前 向 GPS 多 路 径 模型 20。 


P, = PSGiW; (5) 


2 


Per 


XSW, 


(6) 
P 代表 电场 能 量 , G 代表 天 线 增益 , W 代 


表 Woodward 模糊 函数 ， 下 标 4 和 分 别 代 


表 直 射 和 反射 分 量 , X 是 地 表 和 天 线 耦合 系数 。 


X* -R' GF exp(iP) m 
X" -R* VG’ exp(i®’) m 
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脚 标 和 工分 别 为 RHCP 和 LHCP 


极 化 ， 中 为 干涉 相位 ，R* WR 分 别 为 同 
极 化 和 交叉 极 化 发 射 系数 ， 是 H 极 化 和 V 
极 化 的 线性 组 合 。 

4 给 出 了 当地 表 冻 融 转 换 时 ， 天 线 高 
度 为 2.0 m， 载 波 频率 为 GPS LI Bj, KA 
TRM29659.00 天 线 模型 , 模拟 分 析 在 上 半 部 
分 为 空气 , 下 半 部 分 为 冻 融 土壤 的 多 径 反 射 
情况 下 ， 多 路 径 信息 (SNR， 相 位 和 伪 距 ) 
的 变化 情况 ,以 及 其 引起 的 相应 的 有 效 反射 
计 高 度 的 变化 。 从 图 4 可 以 看 出 ， 当 土壤 温 
度 发 生 冻 结 一 融化 转换 时 , 会 引起 多 路 径 幅 
度 增加 ， 以 及 有 效 反射 计 高 度 的 增加 。 而 这 
种 冻 融 转换 引起 的 多 路 径 信息 的 变化 , 可 以 
用 于 地 表 冻 融 转换 监测 。 文 章 [28] 中 给 出 了 
实测 GPS(AB33) 和 SNOTEL (ID 958) 站 点 分 


析 结 果 ， 并 探讨 了 相应 影响 因素 及 局 限 性 。 
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图 4. 土壤 温度 (-0.5°C-0.5°C) 对 GPS 多 路 径 信 息 (SNR、 相 位 和 伪 距 ) 影响 。 图 右 侧 列 图 采用 Lomb 


Scargle 谱 分 析 方法 展示 了 相应 结果 
Figure 4. Soil temperature effects on GPS multipath observations (GPS SNR, carrier phase, and code 
pseudorange). Lomb-Scargle spectrogram analysis is presented in the right panel. Soil temperature with —0.5 °C is 
shown in blue, 0.5 °C in pink. 
2 . 2 冻 融 转化 时 DDM 波形 变化 : a9)» SES BRCS 采用 23 部 分 的 随机 粗 
糙 面 模型 ， 建 立 了 适用 于 冻 融 土壤 的 DDM 
当地 表 发 生 冻 融 转 换 时 , 介 电 常数 的 变 ”模型 , 并 模拟 分 析 了 冻 融 转换 前 后 , DDM 波 
化 会 导致 BRCS 的 变化 , 进而 引起 DDM 波 形 的 差异 。 
形 的 变化 。 本 文 在 Z-V 模型 的 基础 上 (1 ( 公 


A? G 0,9. G 0.9, 
nre [| 4004 0) 
(47) A RR; (9) 
xo, (0.4:0,.6:,0) A" (67) S" (Sf )dA 
Z-V 模型 本 质 上 为 双 站 雷达 的 积分 形 间 的 距离 分 别 用 G 和 R 表示 ， 下 标 工 和 及 


Y (T.f) = 


X. I» GNSS 卫星 发 射 能 量 ， 其 波长 为 ”分 别 代表 发 射 机 和 接收 机 ，C 则 是 积分 区 
4 ， 天 线 增益 及 散射 点 到 发 射 机 /接收 机 之 PA 内 的 冻 融 地 表 的 双关 雷达 散射 截面 , 三 
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MBA Sic 函数 分 别 用 人 A 和 S 表示 ， 二 ”从 中 可 以 看 出 当 冻 融 转换 时 ,不 同 极 化 状态 
的 DDM 马蹄 状 波形 存在 着 明显 的 差异 。 极 
化 信息 是 进行 冻 融 监测 的 重要 指标 信息 , 不 
同 极 化 时 ，DDM 波形 的 最 大 值 和 最 小 值 如 
表 2 所 示 。 

表 2. 线 极 化 和 圆 极 化 时 BRCS 的 最 大 值 和 最 小 值 


者 平方 的 积 是 双 站 雷达 的 模糊 度 函 数 。 
当地 表 发 生 冻 融 转换 时 ，DDM 波形 的 
变化 如 图 5 所 示 GE: 未 进行 归 一 化 处 理 ) , 


Tab. 2 Maximum and minimum BRCS differences at linear and circular polarizations 
Pol VV HH VH RR RL 
DDM differences 
Max 8.85e-18 7.28e-18 0 0 0 
Min -1.54e-17 | -1.40e-17 | -4.57e-18 | -8.40e-17 -9.93e-18 


DDM(F\T_VV) DDM(F\T_HH) DDM(F\T_HV) 
12000 12000 
= 10000 T = 10000 
E] E] 5 
& 8000 & & 8000 
o o o 
a a a 
6000 6000 
4000 4000 
0 100 200 0 100 200 0 100 200 
(a) (b) (c) 
Delay, (GPS L1 chips)*10 Delay, (GPS L1 chips)*10 Delay, (GPS L1 chips)*10 
DDM(FYT RR DDM(F\T_RL 
0 
12000 A 
-2 
X 10000 3 
5 -4 
& 8000 -5 
ia -6 
6000 - 
-8 
4000 9 
0 100 200 0 100 200 
(d) (e) m 
Delay, (GPS L1 chips)*10 Delay, (GPS L1 chips)*10 


图 5. 地 表 冻 融 转 换 时 下 线 各 种 极 化 的 DDM 波形 差异 
Fig.5 The DDM differences at linear (HH, VV and HV) and circular (RR and RL) polarization, when the bare soil 
changes from frozen sate and thawn state. 
3 讨论 依据 。 
本 文采 用 的 冻 融 土 介 电 常数 模型 为 经 
过 实际 验证 的 较为 经 典 的 模型 ,理论 模拟 的 


本 文 从 仿真 模型 角度 分 析 了 利用 
GNSS-R/IR 技术 进行 地 表 冻 融 监测 的 理论 
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结果 与 文章 中 的 模拟 完全 一 致 请 约 ; 各 种 极 
化 的 相干 和 非 相 干部 分 的 双 站 散射 特性 , 在 
模型 验证 时 , 采用 的 方法 是 将 模型 中 椭 倾 角 
和 椭 率 角 修改 为 线 极 化 角度 ,与 原 有 模型 进 
行 比 较 , 结果 一 致 , 因此 本 文 用 圆 极 化 / 圆 线 
极 化 模拟 分 析 冻 融 转 换 时 BRCS 差别 ; 另外 
冻 融 GPS 多 路 径 模型 和 DDM 波形 的 验证 
方法 可 参考 文章 [26,19]。 

极 化 是 电磁 波 重要 特性 , 地 表 反 射 信号 
的 极 化 信息 携带 地 表 重 要 信息 , 极 化 比 是 土 
壤 水 分 反 演 和 植被 状态 研究 的 重要 信息 参 
数 230， 地 表 冻 融 过 程 中 ， 不 同 极 化 的 反射 
言 号 的 双 站 散射 特性 和 DDM 波形 信息 的 差 
异性 表明 极 化 / 极 化 比 信息 在 地 表 冻 融 状 态 
的 监测 领域 具有 重要 前 景 。 

GNSS-R/IR 本 质 上 是 双 站 雷达 
磁 波 散射 的 空 


, 由 于 电 
s 间 异 质 性 , 在 各 个 散射 天 顶 角 
和 方位 角 下 ,包含 地 物 的 不 同 信息 ， 有 效 利 


型 ,同样 为 利用 其 他 数字 通讯 卫星 开展 地 表 
冻 融 监测 提供 了 基本 理论 模型 和 依据 。 
土壤 水 分 、 积 雪 和 植被 都 是 复杂 寒 区 冻 
融 地 表 的 重要 地 表 参 数 ,直接 影响 着 相应 多 
路 径 和 DDM 波形 的 变化 ， 本 文中 给 出 了 裸 
土 情 况 下 ， 地 表 冻 融 状 态 对 GNSS-R/IR 接 


收 机 信号 的 影响 ， 复 杂 寒 区 地 表 的 积 雪 、 植 
被 覆盖 等 是 后 续 模 型 中 发 展 的 重点 。 
4 结论 


GNSS-R 和 GNSS-IR 是 介 于 微波 遥感 
和 卫星 导航 技术 的 新 兴 多 元 交叉 对 地 观测 
技术 ， 近 年 来 在 国内 外 引起 了 广泛 关注 。 本 
文 将 该 种 遥感 方式 的 应 用 领域 从 现 有 的 土 
壤 水 分 、 植 被 和 积 雪 厚度 研究 拓展 到 地 表 冻 
融 状态 监测 上 。 理论 研究 表明 当地 表 发 生 冻 
融 转换 时 ， 土 壤 介 电 常数 存在 较 大 差别 ， 进 


用 不 同 角度 (天 顶 角 和 方位 角 ) 信息 是 提高 
地 物 监测 时 后 向 反 演 的 重点 和 难点 问题 。 
GNSS+R 遥感 是 指 除 了 可 以 利用 GNSS 
导航 卫星 群 的 反射 信号 进行 遥感 监测 外 , JF 
可 以 利用 通讯 卫星 或 数字 广播 卫星 的 反射 
信号 遥感 地 物 参数 11， 本 文 模拟 分 析 时 重 
点 考虑 GPSL1 载波 频率 下 的 散射 特性 、 多 
路 径 以 及 DDM 差异 ， 但 同时 本 文 所 提 及 模 
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而 导致 了 各 种 极 化 下 反射 信号 的 相干 分 量 
和 非 相 干 分 量 在 冻 融 前 后 变化 差异 明显 。 其 
中 ， 多 极 化 相干 分 量变 化 ， 为 利用 GNSS 多 
径 信号 开展 地 表 冻 融 监测 的 GNSS-IR 技术 
提供 了 理论 依据 ; 非 相干 分 量 极 化 状态 的 变 
化 ， 为 利用 漫 散 射 信号 的 DDM 波形 开展 地 
表 冻 融 监 测 的 GNSS-R 技术 提供 了 理论 依 
38. 
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